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We consider class of modified f(R) gravities with the effective cosmological constant epoch at 
the early and late universe. Such models pass most of solar system tests as well they satisfy 
to cosmological bounds. Despite their very attractive properties, it is shown that one realistic 

■ class of such models may lead to significant Newton law corrections at large cosmological scales. 

Nevertheless, these corrections are small at solar system as well as at the future universe. Another 

' realistic model with acceptable Newton law regime shows the matter instability. 

£ 
3 



in 



43 



O 



PACS numbers: ll.25.-w, 95.36.+X, 98.80.-k 



I. INTRODUCTION 



Much attention has been paid recently to the study of modified f(R) gravity (for review, see Such theory which 
may be related with string effective action (H) may successfully describe the dark energy epoch 0, 0, S @, 0, S] • It is 
CIV remarkable that even the form of f(R) gravity may be reconstructed from the known universe expansion history Q. 
Hence, this approach suggests the gravitational alternative for dark energy. It may be considered as proposal for new 
gravity theory which could be more exact than usual General Relativity at current/future universe. If it is so, such a 
theory should pass the known solar system tests 0, []| H, [T(| as well as cosmological bounds. Unfortunately, despite 
the significant progress in the construction of more-less acceptable models, the totally satisfactory theory has not 
been yet proposed. 

Recently, f(R) gravity with the early/late- time effective cosmological constant epoch was proposed @, H, ITU H3|. 
The very attractive, simple versions of such theory [TT|, [l2j seem to show quite satisfactory behaviour from the 
cosmological point of view (the models of ref. 0, [9[ are quite complicated). As well they seem to satisfy (most) of 
solar system tests. Nevertheless, some deviations from General Relativity may be expected. Specifically, the model 
[TT[ may show large Newton law corrections at cosmological scales. Nevertheless, for limited range of parameters 
these corrections are small in Solar System. As well they become small at the future universe. On the same time, 
the model of ref.[T2] may lead to matter instability in the proposed range of parameters. This indicates that such 
theories which show remarkably beatiful behaviour as ACDM cosmologies should be extended, perhaps, introducing 
. ' more parameters. 

II. THE NEWTON LAW CORRECTIONS IN F(R) GRAVITY WITH AN EFFECTIVE 

COSMOLOGICAL CONSTANT EPOCH 

The action of general f(R) gravity (for a review, see (T|) is given by 

S=± J d 4 x^d(R + f(R)) . (1) 

Here f(R) is an arbitrary function. The general equation of motion in /(i?)-gravity with matter is given by 

\g^F{R) - R^F'(R) - g^DF'(R) + V»V V F'(R) = ~T {m)llv . (2) 
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Here F(R) = R + f(R) and T( m \^ v is the matter energy-momentum tensor. 

By introducing the auxilliary field A one may rewrite the action ([T]) in the following form: 

S = i f d 4 x^{(l + f'(A)) (R-A)+A + f(A)} . (3) 

From the equation of motion with respect to A, it follows A = R. By using the scale transformation g^ v — > e° g^ v 
with a = — In (1 + f'(A)), we obtain the Einstein frame action pj]: 



S E = ^ d"x^g <R--( ^tttt g pa d p Ad a A 



= J_ y - \g p °d p ud a o - V(a) 



(4) 



V(a) = e°g ( e -) - e 2 <7 ( ff ( e -)) = ^4) " ^ " < 5 > 

Here g(e~ cr ) is given by solving a = — In (1 + f'(A)) = InF'(A) as A = g(e~ a ). After the scale transformation 
9iiv — * e° <? MI y, there appears a coupling of the scalar field a with the matter. For example, if the matter is the scalar 
field <& with mass M, whose action is given by 

s * = \j d ^V^9 (-g^d^dv* - M 2 <s> 2 ) , (6) 

there appears a coupling with a in the Einstein frame: 

S*e = \J d 4 xV=g(-e a g^d^d u ® - M 2 e 2CT $ 2 ) . (7) 
The strength of the coupling is of the gravitational coupling k order. Unless the mass of a, which is defined by 

(8) 



2 _ 1 d 2 V{a) 



G 2 da 2 

is large, there appears the large correction to the Newton law. 

More exactly, in the Einstein frame, matter fields give a source term for the scalar field a like 

J* = c aa p . (9) 

Here p is the energy density (in the Jordan frame). Now we consider the fluctuations from the background of a = oo 
(co is not always a constant): 

a = cr + Sa . (10) 
For simplicity, we consider the limit where the spacetime is almost flat and consider the point like souces 

Jcr = p^e^^S (x - *«) + p^e a(2> ^5 (x - x™) . (11) 

Then by the propagation of 5a, we find the following correlation function 

/ e a«<T(* (1) V (1 W (3) )^ ^ e (a (1 '+a< 2 )) ( T + a ( 1 )a( 2 )G„( a;i: x 2 ) ^ 

Here G a (x\,X2) is the correlation function of a. When the mass of a is small, we have 

K 2 

G ^ X2 ^ = 12 n(xW -x(2))2 ' (13) 



and 

A(« (1 V (, M* !i) )\ ~ P ( Q(1)+Q(2) ) ff0+ T^CTTy 



(14) 
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At the long range where (x^ — x^) 2 is large enough compared with k 2 , we find 

e (.^)J 1+ fm (2)K L 2 +■■■) ■ (15) 

Then there appears the long range force and the strength of the coupling is given by e( a< >+a< ^""a^'a' 2 '*? 2 . If the 
coupling is very small, the correction to the Newton law might be not so small. 

Recently very interesting f(R) model has been proposed by Hu and Sawicki [TTj | . In the model f(R) is given by 

. N m 2 a (R/m 2 ) n , 

f HS (R) = ,„/ '' , 16 

c 2 {R/m z ) +1 

which satisfies the condition 

lim Jhs{R) = const , 

R— >oo 

Kmf HS {R) = , (17) 

The estimation of ref.fllj suggests that R/m 2 is not so small but rather large even in the present universe and 
R/m 2 ~ 41. Then we have 

f H s(R)~-^ + ^(^Y n , (18) 
C2 c 2 \ mr J 

which gives an "effective" cosmological constant — to 2 ci/c 2 . The effective cosmological constant generates the accel- 
erating expansion in the present universe. Then 

H 2 ~ ??**L ~ (70km/s • pc) 2 ~ (lO- 33 eV) 2 . (19) 

C2 

In the intermediate epoch, where the matter density p is larger than the effective cosmological constant, 

P > 1 , (20) 

C2 

there appears the matter dominated phase (such phase may occur for other modified f(R) gravity as well 7, 8]) and 
the universe could expand with deceleration. Hence, above model leads to the effective ACDM cosmology like models 
0.0- 

Some remark is in order. The approximate expression for the Hu-Sawicky model should be taken with great care. 
The reason is that at very small curvatures where the (non-perturbative) function f(R) goes to zero, the approximation 
breaks down (the corresponding function / may become singular). 

Due to the scalar field in ((4]), an extra (fifth) force could manifest itself. It could violate the Newton law. The 
Newton law is well understood and its correction should be very small at least in the present universe. If the mass of 
a is large enough in the present universe, the problem could be avoided. We now investigate the model by assuming 
A/m 2 = R/m 2 3> 1 since R/m 2 ~ 41 even in the present universe. Then one gets 

nci ( ,4\~" _1 m 2 ci (n+l)m 2 ci ( A\ n 
a —9" ' V ( a ) 2 — 2 ' ( 21 ) 

C2 \ m / C 2 c 2 \ m / 

and 

2 _l d 2 V(a) _ 1 (day 3 d 2 adV 1 / da \ ~ 2 d 2 V _ 1 / A AF(A) 1 | m 2 f A N " +2 

m °^2 da 2 ~~2\dAJ lA2A + 2\dAJ dA 2 ~ 2 \ F'(A) ~ (F'(A)) 2 + F"{A) J ~ 2mTi \m 2 / 

(22) 

First we consider the universe at very large scales, where R ~ (10 -33 eV) and therefore R/m 2 ~ 41. If c\ is not so 
small and/or n is not so large, R/m 2 ~ 41, we find should be very small m„ ~ 10~ 33 eV. Therefore, the correction 
to the Newton law is large. Note that a$ ~ in (fT5|l for the model [ll(. Since a\ t % ~ 1, the correction to the Newton 
law could be not so small. 
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Although m a could be very small at large scales since i?o is very small, Ro can be larger near or in the star. Since 
lg~ 6 x 10 32 eV and 1 cm ~ (2x 10~ 5 eV) , the density is about p ~ lg/cm 3 ~5x 10 18 eV 4 inside the earth. This 

cV) 2 and therefore Ro/m 2 ~ 10 28 . Hence, in 



shows that the magnitude of the curvature could be Ro 



K 2 p ' 
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-19 



case n = 2, we find 



10 GeV, which is very large and the correction to the Newton law is very small. 



Even in air, one finds p ~ 10 6 g/cm 3 ~ 10 12 eV 4 , which gives Ro ~ K 2 p ~ (10 25 eV) and Ro/m 2 ~ 10 16 . In case 
n = 2, m a ~ 10 _1 eV, which gives a correlation length (Compton wave length) about lpm. Thus, the correction to 
the Newton law could not be observed on the earth for such a model. 

What happens in the solar system? In the solar system, there could be interstellar gas. Typically, in the interstellar 
gas, there is one proton (or hydrogen atom) per 1cm 3 , which shows p ~ 10 _5 eV 4 , Rq ~ 10~ 61 eV 2 , and therefore 
R /m 2 ~ 10 4 . Then for n = 2, we find m a ~ 10 25 eV, which corresponds to the correlation length of 10 18 m ~ 100 pc. 
Then the correction to the Newton law could be observed. In case n = 8, however, we find m a ~ 10 _13 eV, which 
corresponds to the correlation length of 10 6 m, which is less than the radius of earth (~ 10 7 m). Then the correction 
to the Newton law could not be observed. Hence, some sub-class of above theory may pass known solar sytem tests 
at the scales of the solar system order. 

In ([22]) . the Einstein frame was considered (J4|. However, similar conclusions may be made also in Jordan frame. 
By multipling §2$ with , one obtains 



3QF'(U) - RF'(R) + 2F(R) 



-T . 



(23) 



Here T = . The equation ([2^)1 corresponds to Eq.(39) in [ll|. Now we consider the background where R is a 
constant R — Rq, that is, (anti-)de Sitter space which can be obtained by solving the algebraic equation 



R F'{R ) + 2F{Rq) = 



(24) 



Since e CT = F'(R), one gets 



SR --F^R) 5<7 



(25) 



Consider the fluctuation 



R — Rq + SR , 



(26) 



which leads to 



1 { F'(Rp) \ __K_ T 



One may consider the point source 



3 V^"(^o) 

T = T 5(x) 



(27) 



(28) 



Then the solution of (f27|) is given by 



da 



k 2 Tq 
6F'(R ) 



G{m 2 ,\x\) . 



(29) 



Here 



(30) 



If m 2 < 0, there appears tachyon and there could be some instability. Even if m 2 > 0, when m 2 is small, 5R ^ at 
long ranges, which generates the large correction to the Newton law. In case of [TTj . we find, when R/m 2 ^ 1 as in 
the present universe, 



Rq 

3n(n + l)ci \m 2 



n+2 



(31) 
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Compared this expression (f3Tj) with (j22[) by putting A — Ro, we find to 2 ~ m 2 . Then the correction to the Newton 
law is the same. 

In it is assumed R/m 2 3> 1 but it might be interesting to study the model assuming R/ra 2 <C 1, which may 
correspond to the future universe. When A/m 2 = R/m 2 <C 1, the potential V(o) (0J) is given by 

V(o) - (1 - n) cim 2 (A/m 2 ) n , (32) 

and we find 

tr~ -ln(l- nc x (A/m 2 )"" 1 ) . (33) 

Let us consider the case n > 1 and < n < 1 separately. 
In case n > 1, when A is small, ()33|) can be written as 



o ~n Cl (A/to 2 )" , (34) 

and therefore 

n/(n-l) 

Then 

-l + l/(n-l) 



V(<t) ~ (1 — n) cito 2 ( —— I . (35) 



n-2 to 
2n(n — l)cj \nci 



(36) 



Note m CT > if ci > 0. Eq. (|34p shows that o is small when A/m 2 is small. Then the mass m a becomes large when 
the curvature R ~ A. Therefore the scalar field does not propagate at large ranges and the Newton law could not be 
violated. 

On the other hand, in case n < 1, for small A, we find 

o (n- l)ln— +ln(-nci) . (37) 

m z 

When A is small, o is negative and large. Eq. ([37|) shows 

V(o) ~ (1 - n) Cl m 2 (-nd) n/(n_1) e^^"" 1 ) . (38) 
In order that the potential being real, ci should be negative. Since 

m 2 „ ( _ nci) n/(n-l) e -WCn-l) (39) 

1 — n 

the squared mass m 2 is negative since c\ < 0, which shows that o is tachyon and unstable. Tachyon is inconsistent 
with quantum theory. Classically if we consider the perturbation with respect to o, the perturbation becomes large. 
Since o is related with the curvature by a = — hiF'(A) = — \nF'(R), the instability may indicate the solution where 
by the perturbation, the curvature of the universe could become large. 

Hence, it seems there may be significant correction to Newton law in the f(R) gravity model under consideration 
at cosmological scales. It is remarkable that such correction becomes negligible in the future, at least, for some range 
of parameters. 



III. THE ABSENCE OF MATTER INSTABILITY 

There may exist another type of instability (so-called matter instability) in f(R) gravity jlCfl. The example of the 
model without such instability is given in [4[ (for related discussions of matter instability, see, [13[ ) . Let us show that 
current and related models are free from such instability. The instability might occur when the curvature is rather 
large, as on the planet, compared with the average curvature in the universe R ~ (lO _33 eV) . By multipling Eq.([2]) 
with g^, one obtains 



(40) 
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Here T = T,j\ . We consider a perturbation from the solution of the Einstein gravity: 



K 2 



R = R o = -—T>0 . (41) 

Note that T is negative since \p\ <C p on the earth and T = —p + 3p ~ —p. Then we assume 

R = R a + R 1 , (\Rx\ « |i? |) . (42) 

Now one can get 



° - ™" + f(W W i?0+ 3F(2)(i? )"3^)(i? )"3^)(i? ) 



, ^ (3) (^o) 
^ 2 )(i?o) 



-□i?! + 2 T?M )Jy p R V p Ri + U(R )R 1 , 



u(Ro) = [Whro] ~ J " * T 

F( 1 Hr )F^(R )Rq _ FW(go) 2F(i? )F( 3 )(fl ) _ FW(R )R 
3F(2)(i? ) 2 3F( 2 )(flo) 3F( 2 ) (R ) 2 3F( 2 )(i? ) 2 



(43) 



If f (i?o) is positive, since Di?i ~ —d 2 R\, the perturbation i?i is exponentially large and the system becomes unstable. 
One may regard V p i?o ~ if it is assumed the matter is almost uniform as inside the earth. 
For the model (fT6|h by assuming Rq/ui 2 ^> 1, it follows 



m 2 Cn 



U W ~ -s^TT) W mJ ) < (44) 

which is large and negative if c\ > 0. Hence, there is no instability in the sense of ref.[l(|. When C\ < 0, however, 
there could be an instability. In first ref. of [l3| , a simple condition for the stability in a sense of [l(| was given, that 
is, theory is stable if F"(R ) = f"(Ro) > but unstable if F"(R ) = f"(R ) < 0. Now we have 

F"(fl )~ " (n+ i )m2ci f^V"' 2 . (45) 



Then F"(Rq) ~ —\/U(Rq) > if ci is positive and theory is not stable. 
As one more example satisfing the conditions (jTTJ) , we now consider 



2 / c 2 («/ m 2 ) r 

m ci ' 



f A {R) = ^^~ e C2( * /m2) " +1 J , (46) 

The asymptotic behaviors of (|4"6"|) are identical with the model (j!6[) when i? is large, 

/a(-R) ~ Ihs(R) -> -— , (47) 

and when i? is small 

~ /ifs(fl) -> ~™ 2 ci (i?/m 2 )" . (48) 

Then asymptotic behaviors of the universe does not change and the correction to the Newton law could be large when 
R is large and small when R is small. The instability is also absent, as one can reobtain the results identical with 
(13211441) . 
Another example is 

f B (R) = -/ e-# , (49) 

with a positive constants fo and to 4 . As in the model (j!6p in (Tlj . we may assume R/ra 2 3> 1 from the early universe 
to the present universe. Even in the model (|47|) . we assume R 2 ^> to 4 . Then by expanding Jb{R) with respect to 
TO 4 /i? 2 , we find 

~ 4 
TO 

MB)~ -fo + fo^ ■ (50) 
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By comparing (|50|) with (fl"8]) . we may identify 

o t m2ci r ~4 ™ 6 ci (M . 
n <-> 2 , /o <-> , / m <-> — — • ( 51 ) 

Hence, /o plays the role of the cosmological constant if /o > 

H 2 ~ /„ ~ (70km/s • pc) 2 - (l(T 33 eV) 2 . (52) 

Thus, the accelerated expansion of the present universe could be generated by the effective cosmological constant /o. 
As in (|20|) . in the earlier but not primordial universe, the matter density p is larger than the effective cosmological 
constant /o- Hence, there occurs the matter dominated phase and the universe could have expanded with deceleration. 
The aymptotic behavior when the curvature is large is identical with the model (|16|) . the correction to the Newton 
law could be not so small. 

We now investigate also the case that the curvature is small. Then for the model (j4"T|) , we obtain 

V(a) = - , + - , - 5 e a* , a ~ — e a* . (53) 
V 2/ m 4 4/ "i 8 y ,4 3 



and 

rf 2 y(q) _ /jaV 3 jV^ /daV 2 rf 2 V A 1 " 4 /A3 

da 2 dA 2 dA \dAJ dA 2 ~ 4f6m 12 / 3 ' 1 j 

If /o is positive, m CT = (l/2)(d 2 V/da 2 ) is positive and large when the curvature R = A is small and therefore there is 
no large corrrection to the Newton law. We should note, however, if fo is negative, which corresponds to the model 
in (UHJ) m 2 becomes negative and there could occur an instability. On the other hand, when the curvature is large, 
U(Rq) in (|4"3"|) has the following form: 



3 

18/qto 

which is negative and large and therefore there is no instability. In fact, since 



U(Ro) ~ -T^Ro , (55) 



f'MRo) ~ 6J jS~ > , (56) 

the condition from first ref. of [HI is satisfied. 

Recently another interesting f(R) model was proposed in [l2j], where 

R I 

Fab(R) = R + fAB(R) = — + — In [cosh(ai?) - tanh(6) sinh(ai?)l , (57) 

2 2a 

with positive constants a and b (for first f(R) models with log-terms, see first ref. in [6j). 

Since the correction to the Newton law has been studied in [12j , we now investigate the possible instability for the 
model ([5T)l . Since 

it is positive. Then the condition [l3| seems to be satisfied and therefore theory seems to be consistent. 
When the curvature R is large, one finds 

F (m R I 1 l, 1 '^ I (l + tanh(6))e- 2afl . 4aR) 

F AB {R)~R+-\n + 2a (l-tanh(6)) + ° {& } ' (59) 

Then U(Rq) (14"3")1 has the following form: 

e^Ro (i _ tanh(fe)) / l-tanh(6)\ , . 

U(Rq) — \- -4t4t l + 21n — ^ . 60 

v ; 6a (l + tanh(6)) \ 2 ) y ' 

If 1 + 2 In ((I — tanh(6)) /2) > 0, U(Ro) is very large and negative and therefore there is no instability. In [12], b is 
choosen to be b > 1.2, so 

I + 21n 1 - ta f (1 - 2) =-3.97<0, (61) 

and therefore the matter instability seems to occur. This indicates that such model should be considered in the other 
range of parameters. 
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IV. DISCUSSION 



In the present letter we considered some solar system tests for several modified gravities which satisfy to conditions 
(16). These theories show very realistic cosmological behaviour and may easily lead to ACDM cosmology. It is shown 
that the theory (15) passes known solar system tests as well as cosmological bounds. Signficant Newton law corrections 
appear only beyond the solar system scales as well as for specific values of curvature power which puts some bound 
for such theory. Theory (59) which has an acceptable Newton law regime shows the matter instability in the proposed 
range of the parameters. Thus, the suggested class of models seems to be very realistic and looks like the alternative 
for ACDM. More accurate and detailed check of cosmological bounds for such theories should be done but in any case 
it is expected that some (combination/extension) of such theories may fit with observable cosmological data. 
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